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Abstract 

An overview is given of different approaches to describing the process of deconfine- 
ment in quantum chromodynamics. The analysis of the known approaches demon- 
strates that the detailed picture of how deconfinement really occurs has not yet been 
understood. Therefore, one has to be rather cautious when interpreting experimental 
signals as attributed to deconfinement. 
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1 Introduction 

When talking about deconfinement in quantum chromodynamics, one should distinguish 
two things: deconfinement as a general phenomenon and deconfinement as a concrete pro- 
cess. What is, in general, the phenomenon of deconfinement is well known - this is the 
transformation of hadron matter into quark-gluon plasma occurring with increasing tem- 
perature or density of matter (see reviews [1-5]). But how this process occurs in reality? 
The answer to the question is crucially important for the correct interpretation of experi- 
mental signals, that are attributed to deconfinement, such as the suppression of charmonium 
production (J/ip suppression [5,6]), the enhancement of strangeness production [7,8], the 
enhancement of dilepton production [7,8] or other signals. 

The process of deconfinement can be considered from two points of view, stationary and 
nonstationary. In the stationary picture, one considers the transformation of equilibrium 
infinite hadron matter into quark-gluon plasma, when rising temperature or density. In the 
nonstationary picture, one takes into account the peculiarities related to heavy ion reactions, 
which involves the consideration of deconfinement in nonequilibrium finite objects [9]. As 
is evident, before going to the complications of the nonstationary picture, it is necessary to 
have a more or less complete understanding of what happens in the stationary picture, which 
the nonstationary one is based on. In what follows, we shall analyze only the stationary 
case and will show that even this has not yet been completely understood. 

2 Numerical Lattice Simulations 

Lattice simulations of high-temperature QCD provide nonperturbative theoretical insights 
into the phenomenology of the transition from hadronic matter to the quark-gluon plasma. 



One of the basic goals of lattice QCD calculations at finite T is to provide quantitative 
results for the deconfinement transition temperature. Up to now, this goal has been achieved 
only in the pure gauge sector [10]. The value of the critical temperature for the deconfining 
phase transition in a SU (3) pure gauge theory (without quarks) is known with small errors 
of the order of 3%. This temperature T c ps 270 MeV, and deconfinement is a real phase 
transition, either of first order or a rather sharp phase transition of second order [10]. 

Unlike in the pure gluodynamics, the transition temperature for SU(3) chromodynamics 
with finite quark masses is not well defined. Even at vanishing baryon number density, there 
is no yet a satisfactory understanding of the critical behaviour in QCD [10]. In zero-density 
chromodynamics at physical values of quarks, deconfinement is rather a rapid crossover than 
a pure phase transition [10,11]. The crossover temperature can conditionally be defined as 
that corresponding to a maximum of some thermodynamic characteristics [12], such as 
a susceptibility [10,11]. Different lattice calculations for the zero-density QCD give the 
crossover temperatures in the interval T c m (140 — 190) MeV [10,11]. 

Finite baryon density calculations in QCD are affected by the so-called sign problem, 
when the fermion determinant becomes complex for nonzero values of the chemical potential 
and the partition function fails to be positive [10,13-15], which prohibits the use of the 
conventional numerical algorithms. Because numerical simulations of QCD at finite baryon 
density are plagued by the principal technical difficulties, the present understanding of 
deconfinement is not satisfactory. However, the available data show no signals of phase 
transition [13-15]. The present situation is rather pessimistic - it seems that there is no 
reliable hope to get important improvements in the knowledge of QCD at finite density 
from lattice simulations [14]. 



3 Pure Phase Models 

Because of the principal difficulties and uncertainty in the finite-density lattice simulations, 
several phenomenological models of deconfinement have been suggested. The most often 
employed are the pure phase models, when hadron matter and quark-gluon plasma are 
treated as different pure thermodynamic phases. Each phase is supposed to possess its own 
thermodynamic potential. 

It is convenient to work with the grand potential Q = —PV, which is a function of 
temperature T, volume V, and chemical potentials fa of different particles, each particle 
sort being enumerated by the index % — 1, 2, . . .. Each type of particles is characterized by 
a set of quantum numbers, such as the baryon number strangeness Si, and others. The 
related baryon density % and strangeness density rig are 

n B = J2 B i Pi n s = J2 S i Pi ' W 

i i 

where pi is the density of particles of type i. Then the chemical potential fa can be expressed 
as 

fa = n B Bi + /j, s Si , (2) 

with the baryon potential fas and strangeness potential ps- Therefore, the grand potential 
can be considered as a function Q = Q(T, V, fas, fas)- Consequently, the pressure is a 
function P = P(T, p B , p s ). The baryon and strangeness densities (1) can be written as the 
derivatives 

dP dP 
n B = t: — , n s = ^ — • (3) 
dfas dps 



In what follows, we shall consider, for simplicity, the case with fixed fis and will analyze 
the dependence of pressure on temperature and fi B , writing P = P(T,/i B ). Respectively, 
the baryon density n B = n B (T, fi B ). 

In pure-phase models of deconfinement, one divides all particles into two groups, one 
group consisting of hadrons and another group consisting of quarks and gluons. The first 
group is assumed to form the hadron phase and the second one, the quark-gluon phase. 
The corresponding pressures of pure phases, Ph and P p , are calculated in different approx- 
imations [1-5]. As a result, the baryon densities of pure hadron and pure plasma phases 
are different, 

^Bh = ^2 Bih Pih = TT-^- , n Bp = ^2 Pv = 7T^~ ' (^) 

i Ofl B . o^ B 

where the summations include, respectively, either only hadrons (and antihadrons) or only 
quarks, antiquarks, and gluons. The deconfinement temperature is given by the equality 

P h (T c , fi B ) = P P (T C , fi B ) (5) 

yielding a uniquely defined transition line T c = T c (fi B ). However, the transition temperature 
as a function of baryon density is not uniquely defined. This is because there are two 
different baryon densities (4), which also are different at the transition temperature, 

n h = n Bn (T c , hm) , n p = n Bp (T c , ^ B ) . (6) 

These densities are different at T c since at this point the pressures of hadron phase and 
of plasma phase intersect. For the deconfinement transition, one has rih > n p . Treating 
the relations (6) as the equations for the baryon potential, one gets two different potentials 
fi h = ii B (T c ,rih) and /i p = fi B (T c ,n p ). Substituting these into the dependence T c (/x#), one 
obtains two lines 

T h (n h ) = T c (fi h ) , T p (n p ) = T c (/i p ) . (7) 

These two lines environ the region on the plane T c — n B where the hadron and plasma 
phases coexist. This situation is completely analogous to the standard case of first-order 
phase transition, when coexisting phases are treated in different approximations or when the 
pressure as a function of density contains an instability interval [16,17]. In the considered 
case, the baryon potential of the hadron-plasma mixture is 

fx B (T c ,n B ) = x h fx h + x p n P , (8) 

where the phase concentrations are defined by the equations 

x h n h + x p n p = n B , x h + x p = 1 , (9) 

which yield 

n B -n p n h -n B 

x h = , x p = . (10) 

n h -n p n h - n p 

Let us stress that the coexisting phases, characterized by the linear combination (8), are 
both macroscopic. Such a mixture of macroscopic phases is called the Gibbs mixture, in 
order to distinguish it from a mixture of mesoscopic phases [18]. 



In the region of the existence of the Gibbs mixture, when n B changes from n p to rih, 
the transition line is represented by a horisontal line 



T c (n B ) = const (n p <n B < n h ) (11) 
connecting the points Th and T p defined in the relations (7). Along this line, the pressure 

p(T,n B )=P(T,fi B (T,n B )) (12) 
is given by a horisontal line connecting the points Ph and P p , according to equality (5), 

p(T c ,n B ) = const (n p < n B < n h ) . (13) 
Then the compression modulus 

K -\ T ,n B )=n B ^- (14) 
on B 

remains zero on the transition line, 

K^ 1 (T c ,n B )=0 (n p <n B <n h ). (15) 

Hence, the compressibility kt = oo diverges exhibiting instability. Thus, the Gibbs mixture 
is, actually, unstable. 

Constructing the pressure of a pure phase, one invokes phenomenological arguments of 
the mean-field type. The most popular is a kind of a quasiparticle description, when each 
sort of particles is characterized by an effective spectrum u>i(k), the interparticle interac- 
tions being included in the spectrum as mean-field parts. For example, one may employ a 
relativistic spectrum 

Ui(k) = ^k 2 + mf + Ui 

or a semi-relativistic one 



UJi{k) = y/k 2 + m1 + Ui , 



where m 8 is a bare particle mass and IT or Ui are the real parts of self-energy playing 
the role of mean fields. The interparticle interactions are taken sometimes in the excluded 
volume approximation, as it is done in the statistical bootstrap models [19] of hadrons. 
The interactions between hadrons can be modelled by various linear [4,12] or nonlinear [20] 
functions of density. 

The quasiparticle picture is also used for the quark-gluon plasma phase [4,12,21,22] 
giving the description of pure gluodynamics and of zero-density chromodynamics in good 
agreement with numerical lattice calculations for the region of pure deconfined phase. How- 
ever, the pure-phase models cannot correctly describe the whole process of deconfinement, 
always predicting a first-order phase transition, contrary to lattice simulations (see discus- 
sion in 



4 Extrapolating Equations of State 

Instead of invoking phenomenological arguments for describing the whole process of decon- 
finement, one sometimes resorts to the following method. One starts with considering the 
properties of a pure phase far from deconfinement, where the equation of state for this pure 



phase can be found with a reliable certainty, and then one extrapolates the found equation 
to the region close to deconfinement. As is evident, analyzing the properties of a sole pure 
phase, it is impossible to get a correct description of the whole phase transition, but only 
the region of stability can be determined in this way. Nevertheless, determining an insta- 
bility point can give a reasonable estimate for that of deconfinement, and the extrapolated 
equation of state could provide an approximation for the region close to deconfinement. For 
example, examining the stability boundaries for nucleons inside nuclear matter yields [23] 
quite reasonable values for the instability temperature of 200 MeV and for the instability 
of 2 normal densities, which are close to the estimates for the deconfinement temperature 
and density [1-5]. 

The QCD pressure is known to be presentable as an asymptotic expansion in powers 
of the coupling parameter g at high temperature and zero chemical potential, when the 
coupling parameter is small. This expansion is known [24,25] to the order of 0(g 6 \ng), 

P(g) -^rT 4 (ao + a 2 g 2 + a 3 g 3 + a A g 4 + a' 4 g 4 \ng + a 5 g 5 ) , (16) 

with the coefficients 

21 / 5 \ / 1 \ 3/2 

a ° = 1+ 32 iV/ ' a 2 = -0.09499 fl + — A^/J , a 3 = 0.12094 M +- N f J 

a A = 0.04331 (l + ^ Nf^j In (l + ^ A^ + 0.01733 - 0.00763V/ - 0.00088A^- 
-0.01323 (l + A ty) (i _ 1 at^ ln|, a' 4 = 0.08662 (l + ± iV>) , 

/ l \ 1/2 

a 5 = - f 1 + - Nf) (0.12806 + 0.00717^ - 0.00027A^) + 

+0.02527 (l + iiV,) 3 ^!- l Nf ) ln£, 

where Nf is the number of flavours. The dimensional regularization is used here, and the 
renormalization scale \x corresponds to the modified minimal substruction scheme MS. 

This expansion (16) is not convergent, but it is merely asymptotic, being valid only for 
g — > 0. Accepting the high-temperature dependence 

g2{T) ~ (HiV c -2iV / )ln(T/A) (T ~" °° ) ' 

in which N c = 3 is the number of colours and A xs 200 MeV is the QCD scale parameter, 
one sees that the condition g <^ 1 corresponds to the very high temperatures T > 10 3 MeV. 
This is why the form (16) does not agree with the lattice simulations for lower temperatures 
[26]. 

To extrapolate expression (16) to the region of finite g, Pade approximants have been 
used [27,28]. The latter are often employed in the attempt of improving perturbative 
results of field theory [29]. However, the constructed Pade approximants exhibit unnatural 
features, containing terms proportional to g both in the numerator and denominator [28]. 
The most important is that the Pade approximants do not converge, but some turn out to 
develop unphysical poles [27,28]. At large g, Pade approximants exhibit chaotic behaviour, 



since Pmn ~ g M ~ N as g — > oo, so that Pmtv — > ±00 if M > N, Pmn — > const for M — N, 
and P M n ^ when M < N. 

Another method of deriving expressions, valid at finite values of the coupling parameter, 
from asymptotic expansions having sense only in the vicinity of zero coupling parameter, 
is based on the Self-Similar Approximation Theory [30-39]. Below, we give a survey of the 
method allowing us to obtain the equation of state in QCD for finite values of the coupling 
parameter and for temperatures in a wide diapason [40]. The approach employs the self- 
similar exponential approximants [38], which, contrary to Pade approximants, contain no 
poles and possess good convergence. 

It is convenient to introduce the dimensionless function 

normalizing pressure P(g) by the Stefan-Boltzmann limit P(0). Then, expansion (16) 
reduces to the set of the approximants 

k 

p k(g) = Em" , 

ra=0 

where k — 0, 1, 2, 3, 4, 5 and the reduced coefficients are 

, _ n _ a 4 + a' 4 lng 
a — 1 , ai = , 04 = , a n = — [n — 2, 6, 0) . 

According to the idea of the optimized perturbation theory [30] , the renormalization scale 
can be treated as a control function defined by the minimal difference condition. For the 
present case, we require that the approximation (18), where \x appears first, be equal to the 
precedent approximation. This leads to the equations 

Pi{g) = PM, 04 = 0. (19) 

Such minimal difference conditions are often employed in theoretical calculations [30,41- 
47]. The meaning of this condition has been explained in the frame of the self-similar 
approximation theory [30-39] as a kind of a fixed-point condition for an approximation 
cascade. Another type of fixed-point conditions is the minimal sensitivity condition [48-55] 
that is also often used in calculations. But the latter condition cannot be directly applied to 
the expansion (18). It is worth noting that the optimized perturbation theory [30] should 
not be confused with the variational minimization of free energy, common in statistical 
mechanics [56,57]. The optimized perturbation theory is a systematic procedure yielding a 
convergent sequence of approximants, while the variational minimization of free energy is a 
one-step procedure giving just a single estimate. In addition, the latter is valid solely for free 
energy or some other thermodynamic potential, since this variation is based on minimizing 
the right-hand side of the Gibbs-Bogolubov inequality, while the optimized perturbation 
theory can be developed for any quantity of interest [30] . 
Condition (19) results in the renormalization scale 

l i = 1 Tg v = l i{T,g) (20) 

as a function of temperature and the coupling parameter, where 

0.01323 (1 + ^^) (l- A At/) hi 7 = 



(18) 



= 0.04331 (l + ^N^j\n(l + ^N f ) 0.01733 - 0.007G3A> - O.OOO88A7 

and 

0.08662 (1 + hNf 



0.01323 (1 + ^iV/j (l-l N f ] 

In particular, for Nf = 6, one has 7 = 0.996964 and v = 5.879155. With the scale (20), the 
form (18) reduces to the expansion 

P(g) ~ 1 + a 2 g 2 + a 3 g 3 + a 5 <? 5 , (21) 

valid for g — > 0. 

The extrapolation of the asymptotic expansion (21) to the region of finite g, by means 
of the self-similar exponential approximants [38,40], leads to the sequence 

P* 2 (g) = exp (c 2 g 2 ) , P* 3 (g) = exp (c 2 # 2 exp(c 3 #)) , 

^5(0) = ex P (c29 2 exp(c 3 g exp(c 5 g 2 ))^ , (22) 

in which the coefficients 

a 2 a 3 a 5 

c 2 = — t 2 , c 3 = — r 3 , c 5 = — r 5 
ao 02 a 3 

are connected with the control functions Tj. The latter are to be defined from fixed-point 
conditions or from the minimization of a cost functional. 

The running QCD coupling a s = a s (fi) satisfies the renormalization-group equation 

» ^ = Pi*.) , (23) 

which, because of the relation g 2 = 47ra s , defines the dependence of g on /i. The renor- 
malization function (3(a) is known for a — > in the four- loop order [58] as the asymptotic 
expansion 

(3(a) ~ b 2 a 2 + b 3 a 3 + b A a A + b 5 a 5 , (24) 

with the coefficients 

If 2 \ , 4 / 19 



6 * = -(i^( 2857 - ir^ + l^ 

65 = - -^-^ (29243 - 6946.3^/ + 405.089NJ + 1.49931iVj) 

In particular, for Nf = 6, we have 

7 52 65 _ 4944.50992 

2^' 3 ~"(4^' 4 ~(4^' 5 ~ (47 



5033 „ 325 



& 2 = -— , b 3 = — — — 777 , 6 4 = — — — , 65 = — = -0.198282 . 



Generally, the signs of the coefficients 6j depend on the number of flavours Nf in the 
following way: 

b 2 < , 63 < , 64 < , 65 < (0 < N f < 5) , 



b 2 < , 63 < , 64 > , 65 < (6 < N f < 8) , 

b 2 < , 63 > , 64 > , 65 < (9 < Nf < 16) , 

6 2 > , 63 > , 64 > , 65 < (17 < N f < 40) , 

6 2 > , 63 > , 64 < , b 5 < (41 < Nf < 00) . 

The qualitative change in the behaviour of (3(a) happens at Nf where the coefficient b 2 
from negative becomes positive [59] . The coefficients b 2 and 63 are renorm-scheme indepen- 
dent, but the higher coefficients 64 and 65 depend on the renorm-scheme employed in their 
calculation. The expansion (24) is obtained [58] within the minimal substruction scheme. 
But since the (3- function does not depend explicitly on fj,, this function is the same in MS 
scheme. 

Defining the reduced function 

13(a) 



m s ^ , (25) 

we find from expansion (24) 

]3(a) ~ 1 + b 3 a + b A a 2 + b 5 a 3 , (26) 

with the reduced coefficients 

b n =^ (n = 3,4,5). 
o 2 

The self-similar exponential approximants extrapolating Eq. (24) to finite a are 

j3l(a) = b 2 a 2 exp(d 3 a) , 131(a) = ^a^exp^aexp^a:)) , 

131(a) = ^a^exp^aexp^aexp^a:))) , (27) 

where the coefficients 

d n = -7^— in (n = 3, 4, 5) 

071-1 

are expressed through the control functions t n , which again have to be defined either from 
fixed-point conditions or from the minimization of a cost functional. 

Substituting the approximants (27) into the renorm-group equation (23), we solve the 
latter obtaining the corresponding approximations for a s (jj). As an initial condition, we 
may take the value 

a s (m z ) = 0.1185 , m z = 91.1882 GeV (28) 

at the Z° boson mass [60]. 

The renorm-group equation (23), with the right-hand side defined by one of the forms 
(27), gives a s = a s (n). The relation g 2 = Ana S) together with the scale (20), leads to the 
equation 

g 2 = A7ra s (fi(T,g)) (29) 

determining the function g = g(T). Substituting the latter in the approximants (22) results 
in the reduced pressure 

p k (T) = P* k (9,T) (30) 

as a function of temperature. 



Calculations show [40] that the behaviour of the reduced pressure p(T) is in reasonable 
agreement with lattice simulations [4]. At the temperature T c m 200 MeV, the pressure 
sharply drops down, when decreasing T, which can be interpreted as confinement. However, 
the details of the confinement-deconfinement process cannot be accurately described by 
such an extrapolation approach based on the consideration of the quark-gluon plasma only. 
In addition, the self-similar exponential approximants (22) or (27) provide an accurate 
extrapolation to the region of finite g or a s , where these parameters are of order one. 
In the vicinity of confinement, when a s fastly grows becoming much more than one, any 
extrapolation procedure would be quantitatively unreliable. 



5 Effective Coupling under Confinement 

The running coupling a s (fj,) as a function of the scale /i is experimentally studied only for 
/i > 2 GeV, where a s < 0.4 [60,61]. In the low-momentum region, the behaviour of the 
effective coupling is poorly known not because of the limited knowledge of higher order 
effects, but because of an essentially different physical phenomenon that enters the game, 
the one that is referred to as confinement. 
For large //, perturbation theory gives 

* ^InWA) ( " " °° ] ' (31) 

where A m 200 MeV is the QCD scale parameter and (3q = — 27r6 2 = 11 — \Nf. The form 
(31) is valid if a s <C 1, that is when 

fi > Aexp > A . (32) 

This implies, as far as A p=* T c ps 200 MeV, that exression (31) is applicable only for p, 3> T c . 
If, nevertheless, one formally considers (31) at lower p, then the coupling (31) diverges at 
p — A. There exist arguments [62] that a s (p) is finite at all p, and satisfies the sum rule 



1 r 2GeV 

- / a s (k) dk ~ 0.38 GeV . 

71 JO 



There are several models [62,63] constructing a s (n) for arbitrary /i. 

The simplest way of a phenomenological construction of finite coupling could be by 
means of the pole-removal trick. The idea of the latter is as follows. Suppose that a 
function f(x) has a pole at x , which implies that in the vicinity of the pole the function 
can be presented as the sum 

f(x) ~ freg(x) + f sin (x) (x -> X ) (33) 

of a regular and singular parts. Let us define a regularized function 

f(x) = f(x) - f sm {x) (34) 

as the function f(x), with the removed singular part. 

In applying this trivial trick to the running coupling, we may postulate that the latter 
is defined by extrapolating the perturbative approximation (31) to all /i by means of the 



regularization procedure (34). The perturbative expression (31) can be identically rewritten 

as 

, . 2nn 
aM ~ /3bln(/i/A)" ' 
with any positive n > 0. Taking into account the asymptotic equality 

11 1 , 

In a; 2 1 — x 

we immediately obtain the regularized coupling 

1 



(x 



2nn 

Po 



+ 



ln(/i/A)« 1 - (/i/A) n 



(35) 



which is finite for any //, including (j, — A and /i = 0, where 

~ / \ 2nn ... 7m . . 

«.(0) = -s- , «.(A) = -=- . (36) 
Po Po 

Similarly, one can construct regularized couplings from perturbative expressions of higher 
orders. But let us note that, as follows from Eq. (35), the pole-removal trick does not 
define the regularized functions in a unique way. 

Another way of regularization, leading to the same result, would be by means of analyt- 
ical continuation. Consider again a function f(x), with x > 0, having one or several poles 
on the positive semiaxis. Assume that f(—x) has no poles. Define the analytic continuation 
f(z) to the complex plane, except the cut along the negative semiaxis, so that 

f(-x±i0)=f(-x±i0) (x>0). (37) 

In the region of analyticity of f(z), the spectral representation 

/>) = - H—dx (38) 
is valid, from where the spectral function is 

J(x) = \ [/(-* + *°) " K-x - M)} ■ ( 39 ) 

In the latter, condition (37) is to be used. 

Applying the analytic continuation method to the perturbative coupling (31), we take 
into account that ln(— x ± iO) = In \x\ ± in and 



In 
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=1= in 


In 2 
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+ 7T 2 



\n(-x ± iO) 
Then the corresponding spectral function is 

Po(m \x\ + n l ) \AJ 

The spectral representation (38) gives the same regularized coupling (35). This way of 
regularization was employed in Refs. [64,65] for the case n = 2. 

The weakest point in any regularization procedure, based on the perturbative expression 
(31), is that the latter is valid only under condition (32), hence a regularized function has 
sense solely for |/x/A| ^> 1. And there is no any reason of extrapolating a s (fi) to the region 
fi/A < 1, where confinement is expected. 



6 Clustering Quark-Hadron Matter 



Any attempt of treating deconfmement from the point of view of pure thermodynamic 
phases contains the following principal contradiction. Hadrons are believed to present bound 
states of quarks and gluons, while quark-gluon plasma represents their unbound states. This 
implies that the system of quarks and gluons possesses, in general, both a discrete spectrum 
corresponding to bound states and a continuous spectrum associated with unbound states. 
When a many-body system possesses an energy spectrum E n , then the distribution of 
particles over the energy levels is described by the Gibbs probability p n ~ e _/3£n , where 
(3T = 1. If the energy spectrum contains both discrete as well as continuous parts, at 
each moment of time there exists a probability for particles to form bound states or to 
pertain to unbound states. In the standard quantum picture, bound and unbound states 
do coexist, with the related probability weights. Hence, hadrons must coexist with quark- 
gluon plasma. Let us stress that this is a direct logical conclusion immediately resulting 
from the treatment of hadrons as describing bound quark-gluon states. If a quantum system 
possesses different parts of spectra, all of them are to be taken into account by calculating 
the corresponding probability weights. It is not correct to separate the whole spectrum 
onto particular sections, prohibiting the existence of some of its parts. It is also incorrect 
to identify separate sections of the quantum spectrum with different thermodynamic phases. 
Hadron states and plasma states are quantum states but not thermodynamic phases. 

This situation is similar to that of electron-ion plasma. Electrons and ions can form 
bound states, i.e. neutral atoms, or unbound states of electrons and charged ions. In the 
system of electrons and ions, under given conditions, there is a fraction of neutral atoms 
and a portion of separated ions and electrons. Changing conditions varies the fractional 
concentrations of neutral and ionized atoms. Ionization in the system of predominantly 
neutral atoms is the direct analog of deconfmement in predominantly hadronic matter. 
Ionization as well as deconfmement can occur, depending on circumstances, either as a 
sharp transition or as a gradual crossover. 

The description of statistical properties of a quantum many-body system possessing 
several qualitatively different quantum states, such as bound and unbound, is not a triv- 
ial task. For describing such systems, Theory of Clustering Matter has been elaborated 
[4,12]. The approach is based on three pivotal concepts: Cluster Representation, Statistical 
Correctness, and Potential Scaling. 

The idea of the quasiparticle cluster representation goes back to the authors who ana- 
lyzed the abundances of chemical elements on Earth by treating each element as a quasi- 
particle characterized by the corresponding atomic weight and the binding energy, with the 
related chemical potentials taking into account the allowed interparticle reactions. Such 
approaches are reviewed in Refs. [66-68]. The same idea was applied to considering nu- 
clear multifragmentation [69]. A more accurate mathematical formulation for the problem 
of constructing the quasiparticle representation for composite particles was initiated by 
Weinberg [70-72]. Such a representation could be unambiguously defined provided that a 
transformation from the state space of elementary particles to that of the system containing 
composite particles, together with unbound elementary particles, would be given [73-75]. 
For this purpose, different Boson realizations of Lie algebras [76] were employed [77-81]. 
The most general approach, based on the Tani transformation [82], has been developed 
by Girardeau [83-86] who coined the term Fock-Tani representation. This was applied to 
various systems containing bound clusters, including the quark-hadron matter [87]. 

The basic point in the quasiparticle cluster picture is as follows. Consider a many- 



body system, with the total space of quantum states being a Fock space T . Let the 
algebra of observables, A, be defined on T. Assume that the particles of the system 
can form several types of bound states, e.g. corresponding to different hadron clusters. 
Enumerate all admissible types of bound states by the index i = 2,3, . . ., reserving the 
index i — 1 to unbound states. Each kind of bound clusters can be individualized by a 
set of characteristic parameters, such as the compositeness number Zi showing the number 
of elementary particles bound into a cluster, effective mass of the cluster rrii, and a set of 
quantum numbers like spin, isospin, colour, baryon number, strangeness, and so on. And 
let us treat each type of bound clusters as a separate sort of particles, with the associated 
Fock space T, called the ideal cluster space. The direct product 

T = ®iT t {T x = T) (41) 

composes the total cluster space. The formal relation between the Fock space of elementary- 
particle states and the cluster space (41) can be presented by means of a unitary transfor- 
mation U, such that 

jr = fjj: ? jr = fj+jr . ( 42 ) 

Then the cluster algebra of observables is defined as 

A ee U+AU . (43) 

With these definitions, all matrix elements of the algebra A in T are the same as those of A 
in T, since TAT = TAT. Since the representations of A in T and A in T are isomorphic, 
all observables quantities are the same in the standard picture of elementary particles and 
in the quasiparticle picture of a clustering system. 

Let us now delineate the mathematical structure of the Tani transformation. Let the 
field operators of elementary particles, say of quarks, be q(x) defined on the Fock space 
T, with x being a set of spatial variables. Suppose ipi(xi,x 2 , ■ ■ ■ , Xi) is a Schrodinger wave 
function describing a bound state of i elementary particles. The field operator of this bound 
state can be presented as 

^i(x) ee J (pi(xi — x : x 2 — x : . . . , Xi — x) q(xi) q(x 2 ) ■ ■ ■ q(xi) dx\dx 2 ■ ■ -dxi . 

The image of the bound state in the ideal cluster space T is given by a cluster with the 
field operator ipi(x). By definition, ipi(x) = q(x). The Tani transformation is described by 
the unitary operator 

U = exp (| F) , P = Hj - *kx)A(x)} dx . (44) 

In the cluster representation, constructed on the cluster space (41), one defines the 
statistical state < A > for the algebra of observables (43). The density of the i- type 
clusters is 

1 f 

Pi = — < Ni > , Ni = J ipt(x)ipi(x) dx . 

The probability for this type of clusters to be formed can be characterized by the weight 



Pi 

Wi = Zi — 



T, Z iPi) > ( 45 ) 



which may be called the cluster probability. The weight (45) satisfies the standard properties 
of probability, being nonnegative, < Wi < 1, and normalized, J2i w i — 1- 
The direct calculation of the cluster Hamiltonian 

H = U + HU (46) 

is a rather complicated problem. Moreover, the actual form of the Hamiltonian (46) is writ- 
ten as an infinite series. Because of this, one usually simplifies the procedure by assuming 
an effective Hamiltonian H e ff, whose construction involves physical reasoning. The latter 
Hamiltonian is often written with an explicit dependence on thermodynamic parameters, 
such as the cluster densities pi and temperature T, so that H e ff = Heff({Pi},T). At this 
point the principle of statistical correctness [4,12] comes into play saying that the general 
form of the cluster Hamiltonian (46) has to be as 

H = H ef f + CV , (47) 

where C is a nonoperator term such that makes the Hamiltonian(47) statistically correct, 
which implies the validity of the equations 

OH dH 

< ^ > = - < w > = - (48) 

From Eqs. (47) and (48) one gets the equations 

dc__± < dth ll> dc _ _ i_ dH^_ 

d Pl ~ V d Pi ' dT V dT > ' 1 } 

defining C = C({pi},T). These conditions garantee the validity of the thermodynamic 
relations 

P = - w = --, £ = T gf -P + P B n B + Psns=y <E>, 

dP 1 dP dP 

s = — = -(£ + P- p B n B - n s n s ) , n B = — = ^ B ^ > n s = a — = S S *P* ' 
oi 1 op B i ops i 

in which e and s are the energy and entropy densities, and 

dP 1 

E = H + PiNi , p i = — = -<N l >. 
i dpi V 

The cluster Hamiltonian (47) contains the terms describing effective interactions be- 
tween different clusters. For defining the corresponding interaction potentials, the principle 
of potential scaling has been formulated [4,12]. According to the latter, the interaction 
potentials from the same class of universality are connected by the scaling relation 



$ij(r) _ $ a6 (r) 



Z{Zj Z a Zfy 



(50) 



This allows the definition of all qualitatively similar interaction potentials through one 
known potential. Another form of scaling (50) could be 

$ij(r) = $q&(r) 

mimj m a mb 



provided that mj ~ Z{. 

The theory of clustering matter has been applied to clustering quark-hadron matter 
[4,12]. The appearance of multiquark clusters in nuclear matter is explained. The possi- 
bility of the dibaryon Bose condensation is advanced. Provisions for nuclear-matter lasers 
are estimated [88]. Thermodynamic characteristics for the SU(3) gluodynamics and zero- 
baryon-density chromodynamics are in good quantitative agreement with lattice simula- 
tions, displaying a first-order transition for pure gluodynamics and a crossover for chro- 
modynamics. Deconfinement at finite baryon density, at conditions typical of heavy-ion 
collisions, is predicted to be a gradual crossover. 

7 Discussion 

The model of clustering quark-hadron matter [4,12] provides, to our mind, the most realistic 
approach to describing deconfinement at finite baryon density. Deconfinement is found to 
be a gradual crossover, but not a sharp transition. However, it would be yet too premature 
to state that all details of the deconfinement process are well understood. We do not 
imply here some technicalities that could always be varied in the frame of the same general 
approach to describing the clustering quark-hadron matter. For instance, one can take 
different interaction potentials or accept different cluster spectra u>i(k). Such technical 
variations do not change the general qualitative picture. But there are more principal 
questions that have not yet been properly addressed: 

(i) The clustering quark-hadron matter has been treated in the mean-field approximation 
[4,12]. This seems to be reasonable especially because deconfinement is not a second-order 
phase transition but rather a crossover. The state where unbound quarks and gluons coexist 
with hadron clusters is shown to be thermodynamically stable with respect to the thermal 
and mechanical stability and with respect to the minimality of the thermodynamic potential 
as compared to those of pure phases corresponding either to quark-gluon plasma or to pure 
hadron states. However, the dynamic stability of the clustering system, which requires 
the positiveness of the collective-excitation spectrum, has not been checked. The latter 
would be suitable to consider as far as in the mean-field approximation the conditions of 
thermodynamic and dynamic stability do not coincide. 

(ii) The finite-size effects in a clustering system have not been analyzed, while this would 
be useful keeping in mind that heavy ions colliding in realistic experiments are always finite. 
The finiteness of a system not only leads to quantitative corrections, as compared to an 
infinite matter, but may sometimes cause the existence of thermodynamic quasi-phases [89] , 
not existing in thermodynamic limit. 

(iii) It is possible that small static bubbles of quark-gluon plasma inside the predomi- 
nantly hadron matter could arise, and vice versa, static hadron bubbles inside quark-gluon 
plasma could exist [90]. Also, static droplets of strange matter, the so-called strangelets 
could be formed. All such possibilities should be considered in the framework of the clus- 
tering matter. 

(iv) Mesoscopic heterophase fluctuations may emerge even in a globally equilibrium 
system [18]. Such dynamically fluctuating germs of quark-gluon plasma or hadron droplets 
are principally different from static bubbles [90] and require a different theoretical approach 
[18]. 

(v) Finally, to be closer to collision experiments, one should consider a nonstationary 



picture, analyzing all different possibilities mentioned above. In the nonequilibrium case, 
several scenarios of deconfinement could be feasible. Then the problem of pattern selection 
would arise, requiring the necessity of defining the probabilistic weights for the admissible 
deconfinement scenarios. 

Summarizing the main material of this review, we come to the following conclusions: 

1. For the correct description of deconfinement, it is necessary to employ the approach 
based on the clustering quark-hadron matter. 

2. The process of deconfinement has not yet been completely understood. 

3. One must be very cautious in trying to interpret observed experimental data as at- 
tributed to deconfinement. 
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